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Abstract. We consider numerically the crossover scaling behavior from the 
directed percolation universality class to the compact directed percolation 
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Our results are compared to those of a recently performed field theoretical 
approach. In particular, the value of the crossover exponent <j> = 2 is confirmed. 
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1. Introduction 

We consider the crossover scaling behavior between two non-equilibrium universality 
classes. Both universality classes describe second order phase transitions into an 
absorbing state where the dynamics of the model is trapped for ever. The first 
class is the famous directed percolation (DP) universality class (see ^ |21 El ED 
for recent reviews). The second class is the so-called compact directed percolation 
(CDP) universality class which is also known as the voter universality class within 
the mathematical literature. Whereas the one-dimensional CDP behavior is well 
understood due to a mapping to random walks, the stochastic process of DP is still 
analytically unsolved. Numerical investigations of this crossover trace back to 
Here, we numerically examine the crossover between both universality classes within 
the one-dimensional Domany-Kinzel (DK) cellular automaton [Hj. This model is well 
suited for the corresponding crossover investigation since its phase diagram exhibits a 
line of DP-like transitions, terminating in a CDP endpoint. In particular, we consider 
the temporal evolution of the order parameter. The behavior of this quantity reflects 
the crossover and will be useful to determine the crossover exponent. Our results 
confirm those of a recently performed field theoretical approach ? . In particular, the 
value of the crossover exponent, which is expected to be exact, is confirmed. 

2. Directed Percolation 

The stochastic process of directed percolation (DP) can be considered as a paradigm 
of non-equilibrium phase transitions into an absorbing state (see [Q 12 [21 0] for recent 
reviews). The process of directed percolation can be represented by the Langevin 
equation [H] (in the Ito sense, see e.g. [Hj) 

\- l d tPa .{r,t) = rp a (r,t) - up^v,t) 2 + V 2 p a (r,i) + V {r,t) (1) 
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which describes the order parameter p a (r,t) 7 i.e., the density of active sites, on a 
mesoscopic scale. A positive order parameter occurs above the transition point (t = 0) 
whereas the absorbing state p a = is approached for negative r, below the transition 
point. Furthermore, 7/ denotes the noise which accounts for fluctuations of the particle 
density p a (r,£). According to the central limit theorem, r/(r, f) is a Gaussian random 
variable with zero mean and whose correlator is given by 

( n(r, t) V (r',t') ) = A -1 T p a (r, t) 5(r - r') S(t - t') . (2) 

Notice that the noise ensures that the system is trapped in the absorbing state 
p a (r, t) = 0. Furthermore, higher order terms such as p a (r, t) 3 , p a (r, t) 4 , . . . (or 
V /9 a (r, t), V 6 p a (r, t), . . .) are irrelevant under renormalization group transformations 
as long as u > 0. Negative values of u give rise to a first order phase transition 
whereas u = is associated with a tricritical point (see ^3 HH E] for field theoretical 
as well as |12II13| for recently performed numerical investigations of tricritical directed 
percolation). 

As usual for second order phase transitions the system obeys certain scaling laws 
close to criticality (see |2] or 0] for complete discussions of the critical scaling behavior 
of absorbing phase transitions). For example, the steady state order parameter scales 
(up to higher orders) as 

p a oc t 13 (3) 

for r > 0. Additionally to the steady state scaling behavior, the dynamical scaling is 
also expressed in terms of power laws 

p a (t) oc t~ a , P a (i) oc t~ s . (4) 

The first law describes the order parameter decay at criticality starting from a fully 
occupied lattice. The second power law describes the temporal evolution of the survival 
probability P a of an initially isolated single seed of activity. This probability is related 
to the probability that a given site belongs to a percolating cluster 

P P c r c oc r' 9 ' , (5) 

for r > 0. Often this quantity is considered as the order parameter of the percolation 
transition, additionally to the steady state density p a . Above the upper critical 
dimension d c DP = 4, the exponents equal their mean field values, e.g. /3 = 1 and 
> I- 

Below the upper critical dimension renormalization group techniques have to be 
applied to determine the critical exponents. In that case path integral formulations 
are more adequate than the Langevin equation approach (see e.g. |14p. Furthermore, 
they provide a deeper understanding of the stochastic process. Stationary correlation 
functions as well as response functions can be determined by calculating path integrals 
with weight exp(— J), where the dynamic functional J describes the considered 
stochastic process. Up to higher irrelevant orders the dynamic functional associated 
with directed percolation is given by [HI 1151 llfil HTj 



J\p^ Pa] = A Jd d rdtp a 



\- l i),l>.. (-<- V> a - ( £p a - U p a )p a 



(G) 



where pjr,t) denotes the response field conjugated to the Langevin noise field |18| . 
The functional J is invariant under the duality transformation (so-called rapidity 
reversal in Reggeon field theory) 

M — 1 Pa(r, t) < — ► -p p a (r, -t) (7) 
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with the redundant variable fj 2 = 2u/T. As usual, the duality transformation defines 
a dual stochastic process that might differ from the original one [T!5] . At criticality, the 
average density of particles of the dual process /?f ual equals asymptotically the survival 
probability ,7, via 

p.(*) ~ M 2 prc*)- (8) 

Here, Pjt) denotes again the probability that a cluster generated by a single seed is 
still active after t time steps. On the other hand, pf ual (£) describes the particle decay 
of the dual process started from a fully occupied lattice. The self-duality of directed 
percolation implies pjt) = /of ual (i) and therefore |2()l I7j 

Pjt) ~ M'PaW. (9) 

The asymptotic equivalence ensures that both quantities have the same exponents |20| , 
including 

a = 5. (10) 

Taking the scaling laws /3 = v^a and j3' = v^8 into account (i/y denotes the exponent 
of the temporal correlations) equation (fTTTfl implies the identity 

P = P (11) 

for the universality class of directed percolation. 

It is worth mentioning that equation (JHJ follows from the rapidity reversal 
symmetry of the dynamical functional J ^ i.e., it is a specific property of the directed 
percolation universality class. Thus, compared to general absorbing phase transitions, 
the number of independent critical exponents for directed percolation is reduced. 
Furthermore, the self-duality is expressed within the field theoretical treatment of the 
minimal Langevin equation of directed percolation. Hence, the important rapidity 
reversal is usually reflected on a coarse grained level only and holds often only 
asymptotically. In other words, it is often masked on a microscopic level, i.e. it does not 
necessarily represent a physical symmetry of microscopic models j? : . Nevertheless, the 
robustness and the ubiquity of the DP universality class is expressed by the conjecture 
of Janssen and Grassberger [HI EI]: short-range interacting systems, exhibiting a 
continuous phase transition into an absorbing state, belong to the DP universality 
class, if they are characterized by a one-component order parameter. Different 
universal scaling behavior is expected to occur in the presence of additional symmetries 
or relevant disorder effects. 



3. Domany-Kinzel automaton and Compact Directed Percolation 

Several lattice models are known which exhibit an absorbing phase transition 
belonging to the directed percolation universality class. Famous examples are the 
contact process [22], the pair-contact process [211], the Ziff-Gulari-Barshad model |24| . 
and the threshold transfer process [2S] (see e.g. j2H] for a detailed scaling analysis 
of various DP-like lattice models). Another well known 1 + 1-dimensional stochastic 
cellular automaton exhibiting directed percolation scaling behavior is the Domany- 
Kinzel (DK) automaton ;0 4 . It is defined on a diagonal square lattice where the lattice 
sites are either empty (n — 0) or occupied by a particle (n — 1). At discrete time steps 
a parallel update procedure is performed according to the following rules (see figure^]) . 
A site at time t is occupied with probability p 2 (pi) if both (only one) backward sites 
(at time t — 1) are occupied. Otherwise the site remains empty. 
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Figure 1. The 1 + 1-dimensional Domany-Kinzel automaton. Occupied sites 
are marked by full circles. A spreading of particles is sketched in the left 
part. The state of the one-dimensional system at time t + 1 is obtained by an 
iteration of the dynamics rules that are embodied in the conditional probabilities 
P(n t +i \ni,t > n-r,t )■ Thus the occupation n t +i of a given site at time t depends 
on the state of the left (1) and right (r) backward neighbors (ni,t,n x ,t) at time t. 
Spontaneous particle creation corresponds to the conjugated field and is usually 
forbidden, i.e., p = 0. 



Notice that the Domany-Kinzel automaton depends on the two parameters p x 
and p 2 . This dependence results in a non-trivial phase diagram which is shown in 
figure|21 Some particular processes are included in the Domany-Kinzel automaton. 
For example, the process of bond directed percolation corresponds to the choice 
Vi = Pi (2 — Pi) whereas site directed percolation is obtained for p 1 = p 2 
Furthermore, forp 2 = the Domany-Kinzel automaton equals to the cellular automata 
rule "18" of Wolfram's classification scheme |2"71l2"5| . 

Instructive insights into the behavior of the system can be gained already from a 
simple mean field treatment (see e.g. [211 GDI El EH)- Within a one-site approximation 
the above presented reaction scheme leads to the differential equation for the density 
of active (occupied) sites (see 0] for a detailed consideration) 

d t p«( Pl ,p 2 ) = 2(2p 1 -l)p. - 2(2p 1 -p a )/£. (12) 

Focusing on the steady state behavior (9 t p a = 0) we find that the order parameter is 
given by 

pJpnPi) = V P«(Pi,p 2 ) = -7T 1 — ■ ( 13 ) 

2^1 - P-2 

The first solution is stable for p ± < 1/2 and corresponds to the inactive phase. The 
active phase is described by the second solution which is stable for p 1 > 1/2. The 
phase diagram of the Domany-Kinzel automaton exhibits two phases separated by 
the critical line p l c = 1/2 (sec figurc|5J). Along this critical line the order parameter 
vanishes linearly (/3 DP = 1) in the active phase 

Pa(px,p 2 ) = -^—Sp + 0(S P 2 ) (14) 

1 -Jh 

with dp — {p t — p c )/p a . Thus, the critical exponent of the order parameter equals the 
mean field value of directed percolation as along as p 2 < 1, i.e., all critical points on 
the transition line (except of the termination point p 2 = 1) belong to the universality 
class of directed percolation. The last statement is also valid for the 1 + 1-dimensional 
Domany-Kinzel automaton although the transition line as well as the values of the 
critical exponents differ from the mean field approximations (see figure^. 

Clearly, the above derived mean field results are not valid for p 2 = 1. In that 
case no particle annihilation takes place within a domain of occupied sites. Thus, 
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creation and annihilation processes are bounded to the domain walls where empty 
and occupied sites are adjacent. This corresponds to the particle-hole symmetry 

n 1 — n (15) 

which is also reflected in the mean field differential equation l)12|l 

d t p& liPa = l) = 2(2^-1)^(1-/0 - 2*pp.(l-p.). (16) 

To be precise, equation i|16|) is invariant under the transformation p a <-> 1 — p a and 
Sp «-> — <5p. Beside the empty lattice, the fully occupied lattice is now another 
absorbing state. The behavior of the domain wall boundaries can be mapped to 
a simple random walk [S] and the domains of particles grow on average for Sp > 
whereas they shrink for Sp < 0. Thus the steady state density p a is zero below p c = l/2 
and p a = 1 above p c . At the critical value the order parameter p a exhibits a jump. 
The associated critical exponent /3 C dp = differs in all dimensions from the directed 
percolation values f3 DP . Since domain branching does not take place the dynamical 
process for p 2 = 1 is often termed compact directed percolation (CDP). The critical 
behavior equals that of the 1 + 1-dimensional voter model ^1 an d it is analytically 
tractable due to the mapping to random walks. Exact results are derived for the 
critical exponents |S1 133] (P = 0, /?' = 1, v\\ = 2, v± = 1), as well as for certain 
finite-size scaling functions [3JD . In particular, the domain growth from a single seed 
exhibits a one-to-one correspondence to a pair of annihilating random walkers. That 
correspondence allows the calculation of the percolation probability [33] 

r o if Pl < 1/2 , 
r>M = {^ if Pl>1/2 , < 17 > 

yielding the percolation exponent /?c DP = 1- hi contrast to directed percolation the 
universality class of compact directed percolation is characterized by the inequality 

/?cop ^ /?cop ■ (18) 

The number of independent critical exponents is therefore greater than for directed 
percolation. In summary, the universality class of compact directed percolation 
is characterized by a continuously vanishing percolation order parameter P per c and 
algebraically diverging correlations lengths 33 , indicating a second order phase 
transition. But due to the misleading discontinuous behavior of the steady state 
order parameter p a the phase transition was often considered as first order. 

Within a field theoretical approach, the CDP process can be described by the 
Langevin equation [7] 

A -1 9 t p s = rp a (l — + V 2 p a + rj . (19) 

In order to ensure that the fully occupied and empty lattice is absorbing the noise 
correlator obeys 

Mr, tHr', t')) = A- X r p.(r, t) [1 - p a (r, t)] S(r - r') S(t - t'). (20) 

Simple dimensional counting shows that the noise is irrelevant for d > 2, i.e., the value 
of the upper critical dimension is <i CjC DP = 2. Associated to this Langevin equation is 
the response functional [7] 



J[p a ,p a ] = d d rdip a 



A- 1 ^ - rp a (l - p a ) - V 2 pa - ^p a p a (l - p a ) 
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Figure 2. The phase diagram of the 1 + 1-dimensional Domany-Kinzel 
automaton. The bold circles indicate the numerically determined transition points 
(the solid line is just to guide the eye). The dashed line corresponds to the simplest 
mean field approximation (see text). In both cases, all transition points along the 
critical lines belong to the universality class of directed percolation except of the 
termination point (pi = 1/2, p 2 = 1). 



The functional is invariant under the transformation 

Pa(r,i) < — > 1 - p a (r, -t) , 
p a (r,t) < — > - p a (r,-t) , 

T < ► - T , 



(22) 



which reflects the particle-hole symmetry and differs obviously from the rapidity- 
reversal symmetry of directed percolation. 

In the following we consider the crossover scaling from CDP to DP. It is instructive 
to express the mean field equation (|12|) in terms of the distances to the critical line 
r = 2Sp and to the termination point re = 2(1 — p 2 ) 

d t p*(pi,p 2 ) = rp a (l-p a ) - Kp\. (23) 

From the CDP point of view (re = 0), the particle-hole symmetry is broken for any 
positive re, i.e., for any vertical distance from the termination point in the phase 
diagram. The crossover from the CDP to the DP behavior can be described in terms 
of scaling forms. For example, the steady state order parameter obeys for any positive 
value of A 

p a (r,re) - A _ft5DP r(TA,/sA*) (24) 

with /3 CDP = and where <j> denotes the crossover exponent. Setting A = kT 1 /^ we 
obtain jO a (r, re) = f(rre -1 /^, 1). Comparing this result with the steady state solution 
(see equation IT3|) 

r/re 



A>(t,k) = 



(25) 



1 + r/re 

we can identify the crossover exponent as well as the crossover scaling function 

x 

4>mw = 1, f MF (x,l) = — . (26) 
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Studying the crossover beyond the mean field level, we have to add to the 
Langevin equation (|19fl an appropriate term which breaks the particle hole-symmetry. 
For example, within the above presented mean field approximation the particle-hole 
symmetry is broken by — Kp\ (see eciuation l23|) . Alternatively, a term of linear order 
p a acts in the same way. The latter case is examined in ref. [7J by Janssen. Performing 
a field theoretical treatment, Janssen derived the crossover exponent 



2 if d<d c = 2 



<P = { d ' (27) 

[ 1 if d>d c = 2. 

Remarkably, these field theoretical results are expected to be exact since the involved 
diagrams can be summed exactly (see [7] for a detailed discussion). In the following 
we consider the 1 + 1-dimensional Domany-Kinzel automaton in the vicinity of the 
termination point and confirm numerically the result 4> d=1 = 2. 



4. Numerical Analysis of the Crossover Behavior 

At the beginning of our analysis, we determine several critical points in order to 
obtain an estimate of the transition line. For this purpose, we consider the dynamical 
properties of the system along various lines p 2 = kp 1 for k ranging from 0.1 up 
to 1.99 (see tableQ]). For each value, we vary p x and determine the critical point. 
Therefor, we use a standard spreading procedure [20| . i.e., we consider the spreading 
of an initial seed and measure the survival probability P.(£). At the critical point 

the survival probability obeys a power-law behavior with the critical exponent 
S = a = 0.159464(6)|21j, see equation (@J . Slightly above or below the critical value, 
the survival probability exhibits deviations from a pure power-law, resulting in a 
significant (left or right) curvature in a double- logarithmic plot (not shown). In this 
way we get a sufficient estimate of the transition line. The corresponding data are 
shown in figure[21 and are listed in tabled 

As well known from scaling theory, the crossover exponent <j> determines, 
additionally to the crossover scaling behavior, the shape of the transition line |36j. 
In our case, the phase boundary is expected to obey asymptotically 

r c oc k 1 ^. (28) 

Here, the crossover parameter k equals again the distance to the termination point 
along the p 2 -direction, i.e., ft oc (1 — p 2 ). On the other hand, r c is given by 

r c oc p 1:C (k) -Pi, (k = 0) = p ljC (p 2 )-i. (29) 

Thus, the transition line behaves as 

p 1:0 (n) = Pl , c (p 2 ) = - + const (1 - p,) 1 ^ . (30) 

Therefore, we plot in figureOHthe crossover parameter 1 — p 2 as a function of p 1(C — 1/2 
in a double logarithmic plot. Approaching the termination point, the transition line 
agrees well with equation (|28|l and confirms the field theoretical result (j> = 2/d 7 , 
i.e., (j> = 2 for d = 1. But as can be seen, equation describes only the asymptotic 
behavior close to the termination point, i.e., for k <C 0.1. This reflects the notorious 
difficulties to determine the crossover exponent <f> from the phase boundary since it is 
in general not clear whether the asymptotic scaling regime is actually reached. 
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Figure 3. The transition points (bold circles) of the 1 + 1-dimensional Domany- 
Kinzel automaton according to equation il-'tOI . The dashed line indicates the 
analytically expected behavior with the crossover exponent <j> = 2. 



As already mentioned above, all critical points along the phase transition 
line (except of the termination point) belong to the universality class of directed 
percolation. Thus the scaling behavior at all determined critical points has to obey 
the rapidity reversal and equation © has to be fulfilled. In order to check its validity 
we have measured first the order parameter decay (started from a fully occupied 
lattice) and second the survival probability (of an initially single seed) at the critical 
point. Examples are shown in figureQ]for p 2 = p t (2 — p x ) (bond directed percolation) 
as well as for p 2 — Pi/2. Bond directed percolation is a particular case since the 
rapidity reversal is microscopically obeyed (see e.g. [2])- Therefore, equation© holds 
for all t with fi = 1, i.e., 

P,(t) = pM . (31) 

This can be seen in figure^] where the survival probability P„(f) and the order 
parameter p a {t) collapse onto the same function. 

Except of this particular case, the order parameter decay equals the survival 
probability behavior only asymptotically. Furthermore, the variable p, is non-trivial, 
i.e., /i^ 1. Both effects can be seen in figure0]for p 2 — Pi/2. After a certain transient 
P a (t) and p a (i) are characterized by the same power-law decay, i.e., a = S. Choosing 
the correct value of /z (here p w 1.4) the tails of both functions collapse onto the 
same curve. Worth mentioning, the transient regime increases by approaching the 
termination point. But more important, the factor /i varies along the transition line. 
Moving along the critical line into the direction of the CDP termination point, fi 
decreases. Approaching the termination point, /i tends to zero signalling the violation 
of equation ©. In other words, /i — * reflects the breakdown of the rapidity reversal 
caused by the change of the universality class from DP to CDP. In figure^! H 2 is 
plotted as a function of the distance from the termination point along the transition 
line. Notice that p 2 vanishes linearly as a function of p x . In that way, p can be used 
to parameterize the critical line of the one-dimensional DK automaton alternatively 
to the crossover parameter n. But in the below presented scaling analysis, we prefer 
to use k since the accuracy of the determination of k, i.e., of the critical line, is 
significantly higher than that of [i. 
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Figure 4. The order parameter decay p a (t) (started from a fully occupied 
lattice) and the survival probability P a (t) (started from a single seed) for two 
distinct critical points of the Domany-Kinzel automaton. For bond-DP the dual 
symmetry is fulfilled microscopically leading to p a (t) = P a (t). In general, this 
equivalence is only asymptotically obeyed i.e., p a (i) ~ M 2 P a (i). Here, the curves 
for p 2 = Pi/2 are shown. 



Table 1. Estimates of critical points of the 1 + 1-dimensional Domany-Kinzel 
automaton. Bond directed percolation {pi =pi(2 — pi)) is denoted by bDP. Site 
directed percolation corresponds to k = 1. The data of 1371 are obtained by series 
expansions. In 1 + 1-dimcnsion, this method yields data of significantly higher 
accuracy than usual Monte Carlo methods. 



k 


Pl,o 


P2,c = fcpi.c 







0.811(1) 







0.1 


0.8015(4) 


0.08015 


1.5747(34) 


0.25 


0.7894(3) 


0.19735 


1.4876(32) 


0.5 


0.7668(2) 


0.3834 


1.3957(31) 


0.7 


0.74515(7) 


0.521605 


1.3195(30) 


0.8 


0.73300(10) 


0.5864 


1.2809(34) 


1.0 


0.70548515(20)|gg 


0.70548515 


1.1902(30) 


1.2 


0.67316(11) 


0.807792 


1.1046(26) 


bDP 


0.644700185(5) 1571 


0.873762052 


1 


1.45 


0.62585(9) 


0.9074825 


0.93737(20) 


1.6 


0.594305(15) 


0.950888 


0.82363(22) 


1.67 


0.57870(8) 


0.966429 


0.76165(26) 


1.8 


0.54865(7) 


0.98757 


0.60868(24) 


1.9 


0.52469(6) 


0.996911 


0.44501(30) 


1.95 


0.5124250(15) 


0.99922875 


0.31626(38) 


1.99 


0.5024969(15) 


0.99996903 


0.14060(44) 



Now we consider the crossover scaling behavior using dynamical simulations. 
Various quantities, for example p a (t) as well as P a (t), reflect the crossover and are 
therefore well suited for our analysis. But in contrast to the survival probability 
(^cdp > 0) the order parameter decay p a (t) has the advantage that the rescaling 
affects only one axis since /3 CDP = a CDP = 0. Thus, for the sake of simplicity, we focus 
our attention to the scaling behavior of the order parameter decay for critical points 
along the transition line. In order to describe the temporal evolution we incorporate 
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Figure 5. Moving along the critical line of the phase diagram, the parameter /i 
decreases (see right figure) until it vanishes at the termination point, reflecting 
the violation of the rapidity reversal for compact directed percolation. Thus /i 
parameterizes the critical line of the one-dimensional DK automaton. 



an additional scaling field to the steady state scaling form i|24[) 

p a (r, t, k) ~ X'PcnP f ( T Aj t \-"u,c D p iK yt>), (32) 

Along the transition line (r = 0), the order parameter has to obey the scaling form 

p a (0,t,/c) - f(0,l, K t^ 2 ) (33) 

where we made used of the analytical results /3 CDP = 0, ^ii^dp — 2 and wlierG wg hcivc 
set t\~ v w^ = 1. 

The inset of figure^] displays the order parameter decay close to the termination 
point (k = 1.8, 1.9, 1.95, 1.99 see Table^). Starting from the fully occupied lattice the 
order parameter remains nearly constant (a CDP = 0) before a crossover to the directed 
percolation power-law behavior takes place (a DP « 0.159). The first transient regime 
increases by approaching the termination point. Using the field theoretical result 
= 20 these different curves have to collapse onto a single curve by plotting p a as a 
function of nt. The obtained data collapse (see figure|f3) confirms the field theoretical 
result. 



5. Summary 

We have investigated the crossover scaling behavior from directed percolation to 
compact directed percolation within the Domany-Kinzel automaton. Our results 
confirm the field theoretical results, in particular the crossover exponent <fi = 2/d. 
Furthermore, we have analysed the rapidity reversal symmetry along the transition 
line. To our knowledge, the important field theoretical parameter /i is numerically 
determined for the first time. Approaching the termination point this parameter 
vanishes, reflecting the breakdown of the rapidity reversal. 

For the sake of completeness we mention that the crossover from DP to CDP 
along the transition line is also reflected in the morphology of the spreading clusters. 
As pointed out in [2], the varying morphology corresponds to the change of short-range 
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Figure 6. The order parameter decay along the critical line of the one- 
dimensional Domany-Kinzel (DK) automaton. Unsealed data is shown in the 
inset (for k = 1.8, 1.9, 1.95, 1.99, from left to right). A data collapse is obtained 
by plotting the order parameter p a as a function of the rescaled time tn oc t(l— p%). 
Both asymptotic scaling regimes (ctcDP = and a D p ~ 0.159) are recovered. 



correlations, leaving the long-range correlations (i.e., the scaling behavior) unchanged 
unless the termination point is reached. 

I would like to thank R. D. Wilhnann for usefull discussions. 



References 

[1] J. Marro and R. Dickman, Nonequilibrium phase transitions in lattice models (Cambridge 
University Press, Cambridge, 1999). 
H. Hinrichsen, Adv. Phys. 49, 815 (2000). 

G. Odor, Rev. Mod. Phys. 76, 663 (2004). 
S. Liibeck, Int. J. Mod. Phys. B 18, 3977 (2004). 

J.F.F. Mendes, R. Dickman, and H. Herrmann, Phys. Rev. E 54, 3071 (1996). 
E. Domany and W. Kinzel, Phys. Rev. Lett. 53, 311 (1984). 

H. -K. Janssen, J. Phys.: Cond. Mat. 17, S1973 (2005). 
H.-K. Janssen, Z. Phys. B 42, 151 (1981). 

N. G. vanKampen, Stochastic processes in physics and chemistry (North Holland, Amsterdam, 
1992). 

T. Ohtsuki and T. Keyes, Phys. Rev. A 35, 2697 (1987). 
T. Ohtsuki and T. Keyes, Phys. Rev. A 36, 4434 (1987). 
S. Liibeck, J. Stat. Phys. 123, 193 (2006). 
P. Grassberger, J. Stat. Mech. , P01004 (2006). 

M. Chaichian and A. Demichcv, Path integrals in physics, Volume 1 (Institute of Physics, 

Bristol, 2001). 
H.-K. Janssen, J. Stat. Phys. 103, 801 (2001). 
C.J. DcDominicis, J. Phys. C (France) 37, 247 (1976). 

H. -K. Janssen, Z. Phys. B 23, 377 (1976). 

P. C. Martin, E.D. Siggia, and H. A. Rose, Phys. Rev. A 8, 423 (1973). 
T. M. Liggett, Interacting particle systems (Springer, New York, 1985). 
P. Grassberger and A. dc la Torre, Ann. Phys. (N.Y.) 122, 373 (1979). 
P. Grassberger, Z. Phys. B 47, 365 (1982). 
T.E. Harris, Ann. Prob. 2, 969 (1974). 

I. Jensen, Phys. Rev. Lett. 70, 1465 (1993). 

R.M. Ziff, E. Gulari, and Y. Barshad, Phys. Rev. Lett. 56, 2553 (1986). 
J.F.F. Mendes, R. Dickman, M. Henkel, and M.C. Marques, J. Phys. A 27, 3019 (1994). 
S. Liibeck and R.D. Willmann, Nucl. Phys. B 718, 341 (2005). 



Crossover scaling in the Domany-Kinzel cellular automaton 



[27] S. Wolfram, Rev. Mod. Phys. 55, 601 (1983). 

[28] G.F. Zcbende and T.J. P. Penna, J. Stat. Phys. 74, 1273 (1994). 

[29] W. Kinzcl, Z. Phys. B 58, 229 (1985). 

[30] T. Tome, Physica A 212, 99 (1994). 

[31] H. Rieger, A. Schadschneider, and M. Schreckenberg, J. Phys. A 27, L423 (1994). 

[32] A. P. F. Atman, R. Dickman, and J. G. Moreira, Phys. Rev. E 67, 016107 (2003). 

[33] J.W. Essam, J. Phys. A 22, 4927 (1989). 

[34] M.J. Kearney and R.J. Martin, J. Phys. A: Math. Gen. 36, 6629 (2003). 

[35] I. Jensen, J. Phys. A 32, 5233 (1999). 

[36] E.K. Riedel and F.J. Wegner, Z. Phys. 225, 195 (1969). 

[37] I. Jensen, J. Phys. A 37, 6899 (2004). 



